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ABSTRACT 

The exquisite angular resolution available with Chandra should allow precision measurements of 
faint diffuse emission surrounding bright sources, such as the X-ray scattering halos created by in- 
terstellar dust. However, the ACIS CCDs suffer from pileup when observing bright sources, and this 
creates difficulties when trying to extract the scattered halo near the source. The initial study of the 
X-ray halo around GX13-hI using only the ACIS-I detector done bv lSmith. Edgar fc Shaferl (|2002D 
suffered from a lack of sensitivity within 50" of the source, limiting what conclusions could be drawn. 

To address this problem, observations of GX13-I-1 were obtained with the Chandra HRC-I and 
simultaneously with the RXTE PCA. Combined with the existing ACIS-I data, this allowed measure- 
ments of the X-ray halo between 2-1000". After considering a range of dust models, each assumed 
to be smoothly distributed with or without a dense cloud along the line of sight, the results show 
that there is no evidence in this data for a dense cloud near the source, as suggested by'Xiang et aLl 
(2005A. Finally, although no model leads to formally acce ptable results, the .W eingartncr & Draine] 
( 2001 ) and nearly all of the composite grain models from IZubko. Dwek fc Arendt I (j2004l) give poor 
fits. 

Subject headings: dust, extinction — scattering — X-rays: ISM 



I. INTRODUCTION 

Practically every band of the electromagnetic spectrum 
affects or is affected by interstellar (IS) dust grains. In 
the IR, PAHs emit lines and small grains emit continuum 
radiation; in the UV/optical, small grains both extinct 
and scatter light. In X-rays, large dust grains (> 0.1/im) 
scatter X-rays, cre ating halos around point sources . 
The classic paper bv lMathis. Rumpl fc Nordsieck I (|1977l . 
MRN77) used the observed optical extinction to deter- 
mine the size distribution of dus t grains between 0.005- 
0.25/im. Newer models, such as IWeingartner fc Draine I 
(|2001i . WDOl), have extended the modeling to include 
polycyclic aromatic hydrocarbons (PAHs) to match the 
observed IR emission as well as other constraints o n grain 
abundances. Recently, IZubko. Dwek fc Arendt I (|2004l . 
ZDA04) found that a wide range of dust compositions 
and size distributions could fit the existing data, and 
suggested that new observational constraints from X-ray 
halos are needed to select amongst these models. 

X-ray dust scattering halos are created by the small- 
angle scattering of X-rays as they pass through dust 
grains. When an incoming X-ray interacts with the elec- 
trons in a grain large compared to the X-ray wavelength, 
the resulting Rayleigh scattering adds coherently in the 
forward direct i on le ading to small-angle s catter ing; see 
Ivande Hulst"" (1957) and iMathis fc Leel (fl99l for de- 
tails, and Draine (20031) for a comprehensive review. 
More generally, the scattering problem can be posed as 
that of a wave interacting wi th a sphere, in which case 
the Mie solution applies (e.g. ISmith fc Dwek Ill998f l. In 
either approach, the scattering depends largely on the 
grain size distribution, with lesser dependencies on the 
grain composition and position along the line of sight. 

Observations of X-ray scattered halos have 
just begun to signif ica ntly impact dust models. 
ISmith. Edgar fc Shafer I (|200l . (SES02)) described 



Chandra observations of GX134-1 with the ACIS-I detec- 
tor and showed that dust grains do n ot have large (> 0.8) 
vacuum fractions considered by iMathis fc Whiffenl 
(,1989) . SES02 also found that the extremely large 
grains found by Ulysses in the solar neighborh ood 
(|Landgraf etall [2000: iWitt. Smith fc Dwek I [Ml do 
not seem to be common throughout the Galaxy. Despite 
these successes, SES02 could not distinguish between 
the MRN77 and WDOl models. This was in part due to 
calibration uncertainties as well as inherent limitations 
of the data. Despite Chandra's excellent angular 
resolution, ACIS-I observations of GX13-I-1 could not 
measure the halo within 50", due to massi ve pileup in 
the A CIS-I detectors. Dr aine . (2003) and iXiang et al.l 
(|2005f ) have both noted that this result is therefore 
insensitive to dust near the source, as scattering from 
dust within the last 25% of the distance would lead to 
features primarily within the excluded 50". To address 
this shortcoming, I obtained a short Chandra HRC-I 
observation of GX13-I-1. The multichannel plate design 
of the HRC-I is far less sensitive to large count rates, 
which allows GX13-|-l's radial profile to be measured 
to within 2" of the source, far closer than previously 
possible. 

2. OBSERVATIONS 

GX13-I-1 was observed simultaneously with the Chan- 
dra HRC-I and RXTE Proportional Counter Array 
(PCA) on February 8, 2005 for 9.1 ksec (ObsID 6093) 
and 6.2 ksec (P90173), respectively. CIAO v3.3 software 
was used to process the Chandra data, which showed 
significant background flares in addition to the flux from 
the bright source. Standard processing was used for the 
RXTE data. 

2.1. Selecting Good Events 
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Fig. 1.— GX13+1 lightcurve from HRC-I (red) and RXTE (in 
blue; divided by 4 for clarity). The full HRC-I field is shown in 
black, and the selected "good time intervals" are shown in green. 
When the full HRC-I spikes much above 100 cts/s, the telemetry 
limit affects the source count rate. 

The full-field lightcurve included significant periods 
when the count rate approached the 184 cts s~^ telemetry 
limit. The expected HRC-I background rate for th e full 
field is - 50 cts s'^ (jChandra X-rav Center Il2006f) . De- 
spite the brightness of the source, the telemetry satura- 
tion was in fact primarily due to the particle background. 
After excluding a 2' radius circle around GX13-I-1, the 
average count rate was < 50 cts s~^, but with excursions 
above 100 cts s~^ where telemetry saturation would af- 
fect the data. To eliminate this problem only time pe- 
riods where the total counts in the field {i.e, > 2' from 
GX13-I-1) were < 45 cts s~^ were included. Although 
this reduced the total good time to 3.58 ksec, ~ 200, 000 
counts were detected within 2' of GX13-I-1 for a source 
count rate of 54.8 cts s~^. All of these effects can be 
seen in Figure [l] which shows the total HRC-I count 
rate, the HRC-I count rate within 2' of GX13-I-1, and 
the RXTE PCA lightcurve (between 2-9 keV) during the 
observation. Within 2" of GX13-I-1 the count rate during 
these 3.58 ksec was 42.7 cts s~^. According to §4.2.3.1 
of the Chandra Proposer's Observatory Guide, the encir- 
cled energy within 2" is w 90%, rising to « 95% within 
10". Based only on these values, it appears that ~ 78% 
of the total counts are "on-axis" while 22% are scattered 
by a combination of interstellar dust and the Chandra 
mirrors. 

The extremely high flux from the source combined with 
the desire to get the highest possible spatial resolution 
required an unusual instrument configuration. In col- 
laboration with the CXC Operations team, the HRC-I 
detector was positioned so that the source would appear 
in one corner of the HRC-I, while still being on-axis to 
the HRMA. This offset retained Chandra's spatial res- 
olution but ensured the source was far away from the 
normal aimpoint. Figure [2[Left] shows the full field of 
the HRC-I, with GX13-f 1 at one corner. 

In Figure [2[Right], a "jet" extending to the NE and 
containing ~ 1000 counts ca n be seen. This j et is a 
well-known detector artifact (|Murrav et al.l [20001 ) which 
is normally removed by the standa rd processing t o a level 
of < 0.1% of the total source flux (jMurrav 112000 ). In the 
case of GX13+1, this jet is ~ 0.5% of the apparent source 
count rate. The most likely cause is the high source count 



rate interfering with the the on-board electronic event 
processing (Dr. Michael Juda, private communication). 
Although the jet could be eliminated with aggressive fil- 
tering, this would also invalidate the standard calibra- 
tion. It was therefore decided to simply ignore all events 
from the "jet" -side of the source, as shown by the box 
region in Figure [^Left] . 

2.2. Extracting the Spectrum and Flux 

The surface brightness of the X-ray halo must be nor- 
malized by the source flux to make absolute measure- 
ments of the dust column density. GX13+1 is observed 
almost constantly since it is a RXTE All-Sky Monitor 
(ASM) source with a average rate of 20-30 cts/s. How- 
ever, as the RXTE ASM has little spectral sensitivity 
and the HRC-I has no effective energy resolution, simul- 
taneous RXTE PCA observations of GX13-I-1 were taken 
to obtain a useful spectrum of the source. Although the 
RXTE PCA itself has only moderate resolution and lit- 
tle sensitivity below 2 keV, GX13-|-l's spectrum is dom- 
inated by 2-4 keV photons (SES02). The PCA spec- 
trum is shown in Figure [3(Left] as fit with a simple 
model consisting of an abs orbed multi- c olor d isk model 
plus a blackbody, following fUeda et aiTl (l2004l). The col- 
umn density was fixed at the value found by Ueda et al. I 
11004) from the Chandra HETG, h = 3.2 x lO^^cm'^, 
since this result is far more accurate than one obtained 
from the PCA. The best-fit inner temperature of the mul- 
ticolor disk and the blackbody were 1.73 keV and 3.52 
keV, with absorbed 1-10 keV fluxes of 7.43 x 10"^ erg 
cm~^s~^ and 2.30 x 10"^*^ erg cm~^s~^, respectively. De- 
spite the large reduced xl (> 50, driven by systematic 
errors), this model is an adequate fit for this work since 
only the total flux and the approximate spectral shape 
are needed to calculate the predicted response of the 
Chandra HRC-I. Nonetheless, it is important to note that 
fluxes measured with the RXTE PCA are systematically 
high by 10-15% in th e 2-10 keV ran ge, compared to other 
X-ray observatories pahoda II2005D . 

To check the expected count rate, we folded this spec- 
trum through the HRC effective area file for on-axis 
Cycle 7 data (hrciD2005-ll-30pimmsN0008.f its), as 
shown in Figure [3[Right]. The total predicted source 
count rate in the HRC-I is 64.6 cts s^^, ^ 18% larger 
than the observed HRC-I count rate of 54.8 cts 
within 2' of GX13-I-1. The discrepancy is primarily due 
to the overestimation from the RXTE PCA calibration, 
with an additional complication due to spatial variation 
in the response of the HRC-I that reduces the effec- 
tive a rea of the detector corne rs relative to the center 
(Donn ellv. Brown fc Hole1 |200^ 

The Chandra PSF, measured as a ratio of the sur- 
face brightness to the source flux, is the background for 
this observation. The RXTE PCA, a non-imaging detec- 
tor, includes both the direct source flux and the scat- 
tered halo photons, which must be removed to avoid 
double-counting. However, as the goal is to measure 
the scattered halo fraction itself, this problem is recur- 
sive. I addressed this by assuming a column density of 
3.2 x 10^^ cm~^ and calculating the total scattered frac- 
tion for the MRN77, WDOl, and ZDA04 BARE-GR-B 
models, weighted by the HRC-I response. The resulting 
halo fraction ranged from 13-26%. This predicted halo 
strength is consistent with the result that 22% of the to- 
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Fig. 2. — [Left] The full HRC-I observation of GX13+1, with the excluded rectangular region marked. [Right] Expanded image of 
GX13+1 on the HRC-I, showing detector "jet" and excluded region. The 2" radius circle shows the near-source region excluded due to 
likely detector non-linearity. 
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Fig. 3. — [Left] The RXTE PCA spectrum of GX13-I-1, fit with an absorbed multicolor disk plus blackbody model. [Right] The best-fit 
model folded through the HRC-I response. The integrated count rate predicted is 64.6 cts/s. 



tal source counts are between 2" — 120". Therefore, for 
purposes of calculating the background PSF, the RXTE 
PCA flux was reduced by 20% to exclude the halo contri- 
bution, with a 7% systematic error. This reduction is in 
addition to the 15% reduction described above. The 7% 
error is likely not the dominant term in the systematic er- 
ror, however. The observation of a very bright source in 
one corner of the HRC-I detector is at the extreme edge 
of the available calibration, and so careful consideration 
of all uncertainties will be required. 

Another concern regarding this observation was that 
a significant short-term change in the source fiux, on 
the order of 12-24 hours, w ould also affect the ha lo in 
a time-delayed manner (e.g. iVaughan et al. |[200l . At 
smaller angles the delay could be even longer. The RXTE 
ASM data was checked for a 10 day period before the 
observation, but no strong or significant variation was 
seen. Although Type I X-ray bursts have been seen from 
GX13-I-1 which s how 3 — 4x the norm al flux, they only 
last ~ 15 seconds (iMatsuba et al.|[l995f ). In this case, no 
bursts were seen in either the HRC-I or PCA lightcurvcs, 
and indeed the halo observation would not be sensitive 
to such a small variation. 



2.3. Point-spread function 

An accurate measurement of the Chandra HRC-I 
point-spread-function (PSF) between 2-100" from the 
source is crucial to this observation. An accurate ray- 
trace model (ChaRT i), of the Chandra HRMA has 
been calibrated for near-source (< 2") photons , SES02 
showed that at large scattering angles this model signif- 
icantly underestimates the PSF, leading to substantial 
problems in the analysis. Therefore, SES02 relied upon 
an ACIS-I observation of Her X-1 as a PSF calibrator, 
but this source is affected by pileup within ^ 10" and 
therefore cannot be used in the 2-10" range. 

As shown in Figure [3[Right] , the spectrum of GX13-)-l 
peaks at ~ 2 keV. The HRC-I's lack of spectral response 
means that spectral differences between any calibration 
source and GX13-I-1 will lead to additional comphca- 
tions. The best possible calibration source would be 
a bright, hard, and lightly-absorbed X-ray source ob- 
served on-axis with the HRC-I. The X-ray binary LJVIC 
X-1 matches these requirements reasonably well, and two 
Chandra observations (ObsID 1200, 1201) of the source 

^ http:/ /asc. harvard.edu/soft/ChaRT/cgi-bin/www-saosac.cgi 
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Fig. 4. — [Top] The radial profile of 3C273's surface brightness, 
divided by the source flux, fit with the sum of a Gaussian plus two 
power laws and a constant. [Bottom] The ratio of the data and 
model, showing small excursions at 1.5" and 10", but generally 
good agreement. 

have been done. However, they were both done early in 
the mission (August 1999) before the HRMA final focus 
was set and are thus unsuitable. Since then, the bright- 
est hard X-ray source with little absorption and a known 
(albeit variable) flux to be observed with the HRC-I is 
3C273 (ObsID 461 on Jan 22, 2000). Figure H shows 
3C273's surface brightness, divided by its source flux, as 
observed with the Chandra HRC-I (excluding the well- 
known jet region). The spectrum was taken from a Chan- 
dra HETG observation done twelve days earlier (ObsID 
459) which is well-fit by an absorbed power-law with T = 
1.67±0.01 andFA'(2-10keV) = (1.08±0.03) x 10-i° ergs 
cm^^s^^. The absorption column was fixed at the Galac- 
tic value, Nh = 1.8 x lO^^cm'^. The predicted HRC-I 
count rate (based on the CXC PIIVIIVIS tool) for this spec- 
trum is 8.7 cts/s, while the actual source count rate was 
26% higher at 11 cts/s. As the source is variable, this 
was taken as showing little change and the flux was sim- 
ply assumed to have increased by 26% during the HRC-I 
observation. A significant but undetected change in the 
spectral shape could also cause this change in observed 
flux and might affect our results. However, BeppoSAX 
observations of 3C273 over a period of 10 days showed 
only small changes in the 2-10 keV spectral shape, with 
r values ranging fr om 1.56 ±0.02 to 1.64 ±0.02 while the 
flux varied by 15% lHaardt et al. I (|1998f) . The potential 
systematic error caused by uncertainty in the true spec- 
trum of 3C273 during its observation with the HRC-I is 
thus smaller than the error due to unavoidable differences 
between the spectra of 3C273 and GX13-I-1. 

The core of the PSF of 3C273 was fit with a Gaus- 
sian term centered at with FWHIVI of 1. 007+0:005" and 
amplitude 1676 ± 42 arcmin^^. In addition, the best-fit 
model included two power-law terms with Fi = 4.06 ± 
0.05, Fa = 2.40+°:g2 and amplitudes Ai = (3.4+J;^) x 
10''^, A2 = (2.12lg J^) x 10-* arcmin-2 at 100". The 
particle and sky background was fit with a constant, 
(2.60 ± 0.02) X 10-3 arcmin-2. As Figure H shows the fit 
is quite good over a large range of surface brightnesses, 
with the somewhat large reduced xt — 2.6 likely due to 
the extreme precision of the measurement compared to 
the relatively simple model. 



3. RESULTS 

The HRC-I observations are most useful between 2- 
100", since beyond that radius the ACIS-I data can 
measure the energy-resolved X-ray halo. Therefore, the 
ACIS-I data were reprocessed (with CIAO 3.3) and rean- 
alyzed following the approach described in SES02 except 
as noted below. Both the HRC-I and ACIS-I results were 
used in the final analysis. We note that in reprocessing 
the ACIS data, the source flux measurement, done via 
the CCD transfer "streak" , was redone with a better cal- 
ibration and improved handling of the background sub- 
traction which resulted in an overall ~ 15% decrease in 
the measured source flux. The calibration changes in- 
clude a spatially-varying modification of order ±5% in 
the quantum efficiency uniformity in CALDB 2.28, and 
an energy-dependent increase of up to 16% in the overall 
effective area which was added in CALDB v3.2.1. 

The data were fit using the CIAO fitting engine Sherpa 
using scattering models ba .sed on the exact Ra yleigh- 
Gans (RG) approximation (jSmith fc Dweklll998[ ). This 
model assumes the grains are spherical but uses the 
energy-dependent optical constants rather than the 
Drude a pproximation wh e n ca lculating the scattering ef- 
ficiency. ISmith fc Dwek I ([1998;) noted that the full Mie 
treatment is necessary for X-rays < 2 keV, since the RG 
overestimates the total scattering at low energies. The 
HRC-I is sensitive to X-rays between 0.08-10 keV, with 
peak efficiency between 0.7-2.0 keV. In all cases the halo 
model for the HRC-I was calculated using an average 
value weighted by the spectrum of GX13-I-1 and effi- 
ciency of the HRC-I. However, since GX13-|-l's spectrum 
as observed by the HRC-I (see Figure [3[ Right] ) is domi- 
nated by photons with E > 2 keV, the use of the simpler 
RG treatment is justified. 

The initial analysis assumed the dust was smoothly 
distributed along the line of sight. Unlike SES02, where 
the predicted PSF was subtracted from the data, here the 
PSF was incorporated into the fitting directly to allow 
for an explicit inclusion of uncertainty in the PSF. Fits 
to the ACIS-I and HRC-I data included a constant fac- 
tor that allowed for calibration uncertainty in the overall 
PSF, caused primarily by systematic errors in the source 
fiux. For both the ACIS-I and HRC-I data this multiplier 
was allowed to vary by up to 10%. In many cases the fit 
pushed the multiplier to an extremum of the range, show- 
ing that systematic uncertainties remain in the data, al- 
though it is not clear what component dominates them. 
In SES02, systematic errors in the PSF manifested as 
energy-dependent column density fits, since to first or- 
der an error in the PSF could be adjusted by changing 
the overall halo scattering. As the total halo intensity 
is inversely proportional to energy, this effect is often 
a linear dependence of the best-fit Nh on energy. To 
check for this, I allowed the value of Nh to vary inde- 
pendently in the HRC-I and each energy band of the 
ACIS-I data. I used the F-test to determine that in only 
one case (ZDA04 BARE-GR-B) were the best-fit ACIS-I 
Nh values better described by a linear energy-dependent 
model than by a constant value. Even in this case, the 
F-tcst significance was only 3.5%, a negligible value given 
the number of different models tried. It seems unlikely, 
therefore, that there is a significant error in the relative 
power in the dust-scattered (halo) and mirror-scattered 
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Table[T]shows the best-fit Nh results for the HRC-I and 
the "average" ACIS-I value fit and the total assum- 
ing smoothly-distributed dust along the line of sight for 
the MRN77, WDOl, and the 15 ZDA04 models. These 
ca n be compare d to th e value of 3.2 x 10^^ cm^^ found 
by lUeda et al.1 (|2004D . The ACIS-I column densities 
are ^ 20% larger than the HRC-I values, although the 
models with the lowest values tend to have the best 
agreement. The most likely cause of this discrepancy is 
cumulative errors in the source flux measurements com- 
bined with calibration differences between the ACIS-I 
and HRC-I detectors. 

Figure [5] shows the profile of the HRC-I data along 
with the ACIS-I data at 2.5 ± 0.1 keV, near the median 
energy of the spectrum as observed by the HRC-I. This 
figure shows the level of agreement between the HRC-I 
and ACIS-I data agree with each other in the overlap re- 
gion (50 — 100"), as well as the large (> 20 x) difference 
in the HRC-I and ACIS-I backgrounds in the 500-1000" 
region. The radial profile shown in Figure [5] is fit assum- 
ing the line of sight (LOS) dust is "smoothly-distributed" 
and has a composition and size distribution described by 
WDOl [Left], and the ZDA04 BARE-GR-S [Right] mod- 
els. Although both models agree with the overall shape 
of the radial profile, the ratio plots show that the WDOl 
model underestimates the ACIS-I data in the 300 — 500" 
range, while the BARE-GR-S model underestimates the 
HRC-I data in the 10 - 50" range. 

The smoothly distributed dust fits show small dis- 
crepancies that might be indications of dusty molecular 
clouds along the line of sight. These would appear as 
"bumps" in the profile whose position and strength de- 
pends upon the relative distance to the cloud and its col- 
umn density. I therefore refit the HRC-I and ACIS-I data 
using a two-component model that included a smoothly- 
distributed component plus a thin cloud with variable 
position and column density; the cloud is treated as a 
sheet with negligible thickness. Models using only a sin- 
gle thin cloud with no smooth component were also con- 
sidered; these gave generally poor fits independent of the 
dust model used, and so were abandoned. This is not 
unexpected since (a) SES02 was unable to fit a single 
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3.3 


COMP-AC-B 


2.9 ±0.2 


2.3 ±0.1 


< 0.0004 


3.6 


COMP-NC-S 


11 4+0-6 


1.3 ±0.3 


< 0.0008 


2.5 


COMP-NC-FG 


^■°-0.2 


1.9 ±0.2 


< 0.0001 


4.3 


COMP-NC-B 




9 =+0.4 


< 0.00001 


6.5 



cloud using only the ACIS-I data and (b) GX13-I-1 is 
reasonably near the Galactic center {{I, b ) = (1 3.5°, 0.1°), 
D = 7 ±1 kpc; iBandvopadhvav et al. I ()1999f )) where a 
sightline dominated by a single cloud would be unusual. 

To reduce fit time, the column density of both the 
smooth component and a cloud was fixed to be the same 
for all datasets, as was the position of the cloud along the 
line of sight. While more realistic than allowing the cloud 
column density to vary as a function of X-ray energy, this 
has the effect of magnifying residual systematic errors. 
As noted previously, the halo strength diminishes with 
energy while the relative PSF strength increases which 
can create a trend in the best-fit column density as a 
function of energy. However, since in only one case out 
of fourteen was such a trend seen previously, it seems un- 
likely that the systematic errors are driving the resulting 
best-fit parameters in the smooth plus cloud model. 

The best-fit parameters for each model are shown in 
Table [51 None of the fits are formally acceptable {xf, 
ranges from 2.1 to 9.9), although some are clearly better 
than others. In cases where the best-fit position is 0, the 
la upper limit is shown. 

4. DISCUSSION 

ISmith. Edgar k Shafer I (|200l) analyzed the ACIS-I 
observations of GX13-|-l's dust scattered halo and found 
that the dust size distribution does not extend to very 
large (> 1/im) grains and that grains do not have a large 
vacuum fraction. However, the ACIS-I data could not 
distinguish between the MRN77 and WDOl models, as 
the two distributions lead to similar scattering profiles at 
large angles. Similarly, the data left open the possibility 
that there migh t be a s ubstantial population of grains 
near the source ([Draine 2003; Xiang et al. 2005). These 
would create a near-source scattered halo that was ob- 
scured by pileup in the ACIS-I detector. The primary 
goal of the HRC-I observation was to remove these un- 
certainties by measuring the halo near the source. This 
would determine which dust model best fit the data, as 
well as detecting (or put limits upon) variation in the 
dust distribution along the line of sight. 

The fit results shown in Tables [T] and [2] contain a few 
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Fig. 5. — [Left] The X-ray halo from smoothly-distributed WDOl 
divided by the source flux. The HRC-I observations arc fit with a red 
HRC-I data to the model is shown below with crosses and the ACIS-I 



surprises. Just as in iSmith. Edgar fc Shafer I ((200l . the 
WDOl model had the smallest column density of any 
of the models when fit with either smoothly distributed 
dust or after adding a cloud. However, the overall re- 
sult was a significantly worse fit than found with either 
the MRN77 or many of the ZDA04 models. As Fig- 
ure [5] shows, the smooth WDOl model fits the HRC data 
well, but underestimates the halo measured by ACIS be- 
tween 150" -400", while the ZDA04 BARE-GR-S model 
underestimates the halo measured by the HRC between 
10 — 50". Examining the other models show that these 
two cases are representative. Adding a single dust cloud 
to the model results in a solution with a cloud 70-90% 
of the distance to the source if the pure smooth model 
underestimates the halo around 30". Conversely, adding 
a cloud component to smooth dust models that underes- 
timate the halo around 300" tend to put the cloud near 
the Sun. 

Although uncertainties remain due to calibration is- 
sues, the overall quality of the fits shown in Table [1] and 
Figure [5] do not support the proposition that a signifi- 
cant cloud of dust is present near G X13-I-1. In particu - 
lar, although it is not identical to the lXiang et al.l ()2005f ) 
model. Figure [5] shows the result from a model similar 
to their fit to the zero-order HETG for GX13+1 using 
WDOl-type dust. This approximation to their model 
puts a cloud with Nh ~ 2.7 x 10^^ cm^^ at a position 
99.5% of the distance to GX13-I-1, along with a smooth 
distribution with Nh — 1-81 x 10^^ cm^^. In either their 
model or my approximation, the cloud near the source 
dominates the variation in the dust distribution. Fig- 
ure [H] shows that although the HRC-I data contain an 
obvious halo, this model is not a good fit and, in fact, re- 
moving the cloud improves the fit. Although their model 
fit can be improved somewhat with a small vertical off- 
set (corresponding to different relative flux calibrations 
of ^ 12%), the model still predicts 10-20% more halo for 
9 < 10" than is observed. Table [2] confirms this con- 
clusion, finding (in the case of WDOl-type dust) a cloud 
near the Sun, not the source, with Nh ~ 5 x 10^^ cm~^. 

ZDA04 described in detail how constraining dust mod- 
els requires combining multiwavelength data from the 
IR to X-rays while simultaneously considering the metal 
abundances in the grains. Due to the nature of op- 
tical/UV extinction and X-ray scattering, few sources 



type grains fit to the radial profile of GX13-|-l's surface brightness 
line; the ACIS-I data (at 2.5 keV) with a blue line. The ratio of the 
data with circles. [Right] Same, for the ZDA04 BARE-GR-S model. 

show strong signatures of dust in all of these wavebands 
(jValencic fc Smitli1l2007D . Nonetheless, it is possible to 
constrain the allowed dust models by comparing the col- 
umn density predicted by the models to that measured 
using other techniques. In the case of GX13-fl, mea- 
surements of the column density ran ge from 2.5 — 4.0 x 
10^^ cm"^ (iCharles fc Navlor|[l99^ . Optical measure- 
ments provide only an upper limit of 2.9 x 10^^ cm~^ 
based on plausible b ut unconfirmed assum ptions about 
the source spectrum ()Garcia et al. I |1992[ ). The total 
H I column density through the Galaxy at the posi- 
tion o f GX13-hl is 1.8 x 10^2 cm-2 pickev fc Lockman I 
IT99O), but this is misses the contribution from molec- 
ular H2 that is likely to be substantial in the Galac- 
tic plane. The HETG ob servation of GX13 -I-1 agrees 
(weakly) with these results (lUeda et al. 112004 ). although 
it does not strongly limit it. Only Mg can be directly 
measured {Nj^[g — 1.84^049 x 10"^* cm"^), equivalent to 
Nh= 4.81^ 3 ^ ^^^^ cm~2 assuming solar abundances. 
The 2(7 upper limits for Si and S are equivalent to 
Nh < 4 X lO'^^ cm~2. Of course, LMXBs have shown sig- 
nificant variable internal absorption (Hertz fc GrindlavH 
[1983^ inX -rays, so this spectral measurement sets at best 
an upper limit to the actual interstellar component that 
is responsible for the halo. 

Despite these difficulties in independently measuring 
the total LOS dust column density, we can reasonably 
justify excluding the value of Nh = (1.11 ± 0.06) x 
10^3 cm-2 found in Table ffl for the ZDA04 COMP-NC- 
S model fit to the HRC-I data. However, this was only 
model from ZDA04 that used composite grains without 
bare carbon grains that had a plausible value of xt- 
though more data are needed, this class of models, along 
with the group of "composite grains with bare amor- 
phous carbon" models are clearly suspect since they do 
not generate an X-ray halo similar to these observations. 
In fact, the only smoothly-distributed ZDA04 composite 
grains model that fit with < 3 had graphitic carbon 
and B star abundances (COMP-GR-B). After adding a 
dust cloud to the model, the COMP-GR-B model fit with 
x1 = 2.2 while the next best fit (excluding the unrealistic 
COMP-NC-S model) was COMP-GR-FG with xl = 2.6, 
a significantly worse fit. 

5. CONCLUSIONS 
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Fig. 6. — [Top] Ra dial profile of GX13 +1 from the HRC-I (upper points) with PSF calibration data from 3C273 (lower points). A close 
approximation to the lXiang et al.l 1)20051 ) model for GX13+1 using the WE)01 model is shown in blue, with the PSF fit to the 3C273 data 
is shown in red. [Bottom] Ratio of GX13+1 data to lXiang et al.l (120051 ) model. 



The principal results from this analysis are: 

1. Although challenging, HRC-I observations can be 
used to recover the near-source region excluded by 
pileup in the ACIS-I detector. The lack of energy 
resolution can be finessed if another measurement 
of the source spectrum is available. 

2. Fitting the source profile and background PSF in- 
dependently improves overall results, since calibra- 
tion uncertainties in the flux from the source and 
background objects can then be included explicitly. 

3. There is no strong signature of a dust cloud at 
the source i i i the radial profile, as suggested by 
IXiang et all (|2005f) . although some models include 
a cloud ~ 90% of the distance to the source. 

4. In agreement with SES02, the WDOl model under- 
estimates the total column density to the source, 
and again leads to poor fits, although not so bad 
as to be excluded given the calibration uncertain- 
ties. 

5. Some models from the ZDA04 paper, if not con- 
clusively excluded, are at least implausible. In 
general, the ZDA04 models with composite grains 
(excepting the graphitic carbon model with B star 
abundances) gave poor fits, while the bare carbon 
and silicate grain models tended to fit well. 

It should be noted that the relatively good fits found 
using the simple smoothly-distributed dust model are 



somewhat surprising, since X-ray halos probe both the 
largest grains whose size and composition are the least 
constrained from observations in other wavelengths. Ad- 
ditionally, X-ray halos are primarily observed through 
highly- absorbed lines of sight. These probe dust in dense 
molecular clouds that cannot be observed in the optical 
or UV due to the extremely large extinction. Finally, 
all of these models assume spherical grains, although re- 
cently some calculations have been done on aspherical 
grains (jPraine fc AUaf-Akbaril l2006l ) that show the ef- 
fects are small (< 20%) for the WDOl model at ~ 2keV. 
Despite these potential problems, a number of existing 
grain models agree quite well with the observations, sug- 
gesting grains in dense clouds (with the exception of the 
densest regions that take up very little volume and may 
be optically-thick to X-rays) are not too dissimilar from 
grains in less dense regions. 
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